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In this paper the implications of a recently proposed phenomenological model of cosmology, the 
Asymptotic Cosmological Model (ACM), on the behavior of scalar perturbations are studied. Firstly 
we discuss new fits of the ACM at the homogeneous level, including fits to the Type la Supernovae 
UNION dataset, first CMB peak of WMAP5 and BAOs. The linearized equations of scalar pertur- 
bations in the FRW metric are derived. A simple model is used to compute the CMB temperature 
perturbation spectrum. The results are compared with the treatment of perturbations in other 
approaches to the problem of the accelerated expansion of the universe. 
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INTRODUCTION 

According to General Relativity (GR), if the universe 
is filled with the particles of the Standard Model of par- 
ticle physics, gravity should lead to a decelerated ex- 
pansion of the universe. However, in 1998 two indepen- 
dent evidences of present accelerated expansion were pre- 
sented [H [2] and later confirmed by different observations 

ISIIIS]. 

There is no compelling explanation for this cosmic ac- 
celeration, but many intriguing ideas are being explored. 
These ideas can be classified into three main groups: new 
exotic sources of the gravitational field with large nega- 
tive pressure |6j (Dark Energy), modifications of gravity 
at large scales and rejection of the spatial homogeneity 
as a good approximation in the description of the present 
universe |8]. 

Different models (none of them compelling) for the 
source responsible of acceleration have been considered. 
Einstein equations admit a cosmological constant A, 
which can be realized as the stress-energy tensor of empty 
space. This A together with Cold Dark Matter, Standard 
Model particles and General Relativity form the current 
cosmological model, ACDM. However, quantum field the- 
ory predicts a value of A which is 120 orders of magni- 
tude higher than observed. Supersymmetry can lower 
this value 60 orders of magnitude, which is still ridicu- 
lous |9]. In order to solve this paradox, dynamical Dark 
Energy models have been proposed. 

This has also lead to explore the possibility that cos- 
mic acceleration arises from new gravitational physics. 
Again here several alternatives for a modification of the 
Einstein-Hilbert action at large and small curvatures jlO) . 
or even higher dimensional models |1H [12], producing 
an accelerated expansion have been identified. All these 
analyses include an ad hoc restriction to actions involving 
simple functions of the scalar curvature and or the Gauss- 
Bonnet tensor. This discussion is sufficient to establish 
the point that cosmic acceleration can be made compati- 
ble with a standard source for the gravitational field but 
it is convenient to consider a more general framework in 



order to make a systematic analysis of the cosmological 
effects of a modification of general relativity. 

The Asymptotic Cosmological Model (ACM) was pre- 
sented [13] as a strictly phenomenological generalization 
of the Standard Cosmological Model including a past and 
a future epoch of accelerated expansion. It follows from 
the assumptions that GR is not a fundamental theory, 
but only a good approximation when the Hubble rate H 
is between but far away from two fundamental scales 
and H+, which act as bounds on H. A general covari- 
ant metric theory of gravity without spatial curvature is 
assumed. The model is well defined in the homogeneous 
approximation and includes ACDM as a particular case. 

In next section we review the ACM and we provide 
new fits to the Type la Supernovae UNION dataset, first 
acoustic peak of CMB of WMAP5 and BAOs. In the 
third section we derive the linearized equations of the 
scalar perturbations of the metric, following from general 
covariance and a single new assumption on the perturba- 
tions. In the fourth section we will consider how to solve 
the system of equations for adiabatic perturbations in a 
given fluid. In the fifth section we derive the CMB spec- 
trum in the ACM. In the sixth section we will compare 
the treatment of the scalar perturbations in the ACM 
with other models which are equivalent in the homoge- 
neous approximation. The last section is devoted to the 
summary and conclusions. 

THE ASYMPTOTIC COSMOLOGICAL MODEL: 
HOMOGENEOUS BACKGROUND 

The Asymptotic Cosmological Model (ACM) was in- 
troduced in Ref. [13]. In this model the universe is filled 
with photons and neutrinos (massless particles), baryons 
(electrically charged massive particles) and Dark Matter 
(electrically neutral massive particles), but General Rel- 
ativity (GR) is only a good approximation to the gravi- 
tational interaction in a certain range of the Hubble rate 
H, between but far from its two bounds, H- and 
General Covariance and the absence of spatial curvature 
are assumed. 
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The gravitational part of the action might include 
derivatives of the metric of arbitrarily high order, and 
therefore arbitrarily high derivatives of the scale factor 
should appear in the Friedman Equations. We should 
start by considering a generalized first Friedman equa- 
tion 



87rGp(i) = 3/(i/(t),ij(i),...,i/(")(i),-) 



(1) 



and the corresponding equation for the pressure will be 
derived using the continuity equation 



(2) 



However, as the resulting differential equations should 
be solved and only one of its solutions deserves interest 
(the one describing the evolution of the universe), we 
can use the one to one correspondence between time t 
and Hubble parameter H (assuming < 0) to write 
the modified First Friedman Equation as a bijective map 
linking the total energy density p with H in our universe 



(3) 



The use of the continuity equation ([2| enables us to write 
the modified Friedman Equation for the pressure evalu- 
ated at the solution corresponding to the cosmic evolu- 
tion 



SttGp = 3g{H) + g'{H)H/H . 



(4) 



The evolution of the universe will be determined by 
the concrete form of the function g{H), which we assume 
to be smooth. However, the most significant features of 
the evolution at a given period can be described by some 
simple approximation to g{H) and the matter content. 
Those are a pole at H = of order a+ and a zero 
of order q;_ at iJ = The energy density of the 

universe has a contribution from both massless particles 
(radiation) and massive ones (matter). 

For the sake of simplicity, the history of the universe 
can be then divided into three periods. In the first pe- 
riod, H ^ H_ so the effect of the lower bound can be 
neglected and the universe is radiation dominated. The 
universe undergoes (and exits from) an accelerated ex- 
pansion which we call inflation. The simplest parameter- 
ization is 



Hi 



(5) 



In the second period we cannot neglect neither the ef- 
fect of radiation nor the effect of nonrelativistic matter, 
but ^ H ^ H„ . Then GR offers a good description 
of the gravitational interaction, g{H) w in this re- 
gion, and the universe performs a decelerated expansion. 



In the third period, i7+ ^ H so the effect of the up- 
per bound can be neglected, and the universe is matter 
dominated. This period corresponds to the present time 
in which the universe also undergoes an accelerated ex- 
pansion. The simplest choice is 



9iH) ^H' I 



Hi 



(6) 



We will assume that the transitions between these 
three periods are smooth and that the details about these 
transitions are unimportant. 

This model preserves the successes of the Cosmological 
Standard Model, while giving a description of the early 
accelerated expansion (inflation) and of the present one, 
including ACDM as a particular case (a_ — 1). 

Without any knowledge of the evolution of perturba- 
tions in this model, the background evolution at late 
times can be used to fit the Type la Supernovae UNION 
dataset [T3], the first acoustic peak in the Cosmic Mi- 
crowave (CMB) Background p!5], and the Baryon Acous- 
tic Oscillations (BAOs) p[6], via the parameters Hq (the 
present Hubble parameter), H- and a-. 

We use Monte Carlo Markov Chains to explore the 
likelihood of the fit of the supernovae UNION dataset 
|14] to the ACM. The dataset provides the luminosity 
distance 



di{z) = c(l + z) / dz'/H{z') 



(7) 



and redshift of 307 supernovae. A analysis have been 
performed, where has been marginalized over the nui- 
sance parameter Hq using the method described in |17] . 
The resulting parameter space is spanned by the values 
of a_ and H./Hq (FIG. 1). 

The constraints from the CMB data follow from the 
reduced distance to the surface of last scattering at z = 
1089. The reduced distance R is often written as 



R = nll^Ho / dz/H{z) 
Jq 



(8) 



The WMAP-5 year CMB data alone yield Rq = 1.715 ± 
0.021 for a fit assuming a constant equation of state u; 
for the dark energy [21] . We will take this value as a first 
approximation to the fit assuming ACM. We can define 
the corresponding as = [{R — Ro)/'^Ro\'^j ^^'^ ^^'^ 
the confidence regions of the joint constraints (FIG. 2). 

BAO measurements from the SDSS data provide a 
constraint on the distance parameter A{z) at redshift 
z = 0.35, 



A{z) = nl/'HoH{z)-^ 



-1/3^-2/3 



dz'/H{z') 



2/3 



(9) 



Ref. |T6] gives Aq = 0.469 ± 0.17. We can define the cor- 
responding as = [{A{z = 0.35) - AQ)/uAof ■ The 
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FIG. 1: Confidence regions in parameter space of the Asymp- 
totic Cosmological Model (ACM) from the fit of the super- 
novae UNION dataset without priors at la, 2a and 3a. The 
ACDM is inside the la region. The best fit to ACM lies in 
a- = 0.35, H-/Ho = 0.97 (x^ = 310.5) in contrast to the 
best fit to ACDM, which lies in q_ = 1 , H^/Hq = 0.84 
(x' = 311.9). 
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FIG. 2: Confidence regions in parameter space of the Asymp- 
totic Cosmological Model (ACM) from the fit of the super- 
novae UNION dataset (red), of the distance to the surface of 
last scattering from WMAP-5 (blue) and from the joint fit 
(magenta) at la, 2a and 3a. Significantly, the ACDM is still 
inside the la region, unlike in our previous study. This is due 
to the change in the Supernovae dataset. The best fit to ACM 
lies in a_ = 1.50, H-/Ho = 0.77 (x^ = 312.5) in contrast to 
the best fit to ACDM, which lies in q_ = 1, H-/Ho = 0.86 
(X' =313.2). 



confidence regions resulting from adding this constraint 
are shown in FIG. 3. 

In our previous work the use of the supernovae Gold 
dataset |^ and of the WMAP-3 data [H] led us to the 
conclusion that the ACDM was at Scr level in the param- 
eter space of the ACM. The position of the confidence 
regions seems to depend very tightly on the dataset that 



FIG. 3: Confidence regions in parameter space of the Asymp- 
totic Cosmological Model (ACM) from the fit of the super- 
novae UNION dataset (red), of the distance to the surface of 
last scattering from WMAP-5 (blue), of the Baryon Acous- 
tic Oscillations peak (green) and from the joint fit (grey) at 
la, 2a and 3a. Measurements of the BAOs peak do not add 
significant information with their present precision. The best 
fit to ACM lies in q_ = 1.50, H-/Ho = 0.77 (x^ = 312.7) 
in contrast to the best fit to ACDM, which lies in a_ = 1, 
H-/Ho = 0.86 (x^ = 313.7). 

is being used. However, the value of the combination of 
parameters 

(10) 

does not depend much neither on the value of a- or the 
dataset used (FIG. 4). 

Moreover, we can conclude from the figures that BAO's 
do not provide much information in order to constrain the 
confidence regions of the ACM, unlike in other models 
such as ACDM with nonzero spatial curvature. 

SCALAR PERTURBATIONS 

The lack of an action defining the ACM is a serious 
obstacle in the derivation of the equations governing the 
behavior of the perturbations. Given a background be- 
havior described by the ACM, what can be said about 
the evolution of perturbations on top of this background? 
We will find that general covariance together with an ad- 
ditional assumption fixes completely the set of equations 
for the scalar perturbations in the linearized approxima- 
tion and in the region close to 

The key point will be the following. Knowing the exact 
Friedman equations in the homogeneous approximation 
gives us a clue on the form of the equations for the scalar 
perturbations. In particular, we can formally describe 
perturbations over the FRW metric which do not de- 
pend on the spatial coordinates, {0(x, t) = (j){t),ip{'K, t) = 
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FIG. 4: Confidence regions in parameter space of the Asymp- 
totic Cosmological Model (ACM) from the fit of the super- 
novae UNION dataset (red), of the distance to the surface of 
last scattering from WMAP-5 (blue), of the Baryon Acoustic 
Oscillations peak (green) and from the joint fit (grey) at la, 
2a and 3a. The plot shows a- in the horizontal axis and 
flm in the vertical axis. Noticeably, Q,m = 0.26 ± 0.04 almost 
independently of the value of a_. 

ip{t)}. Thus the perturbed metric becomes the FRW 
metric written in a new coordinate frame. With a change 
of coordinates one can derive, starting from the Friedman 
equations, the terms containing only time derivatives of 
the scalar perturbations. Next an assumption on the va- 



SnGSp 
87rG(p(o) +P(o))5j6i 
8TTG{-5p5} + didjU) 



The coefficients anmT--,hnm are adimensional func- 
tions of the Hubble parameter an its time derivatives, 
and can be turned into functions of just the Hubble pa- 
rameter using the bijection explained in the beginning of 
the previous section. We have mentioned we are going 
to be able to determine exactly the terms with only time 
derivatives, that is, the precise form of the coefficients 
o-Ormbom, eom s-iid /om- This shows that there is much 
freedom of choosing a covariant linearized theory of cos- 
mological perturbations, even for a given solution of the 



lidity of the GR description of scalar perturbations when 
H H+, together with the relations, valid for 
any general covariant theory, between terms with time 
derivatives and those involving spatial derivatives, allow 
to derive the evolution of scalar perturbations. 

After this short sketch, we will perform the derivation 
of the equations for the scalar perturbations in detail. 
We can write the metric of spacetime with scalar pertur- 
bations in the Newtonian gauge, 

ds^ = (1 + 20(x, t))dt^ - a2(t)(l - 2V'(x, t))dii^ , (11) 

and the stress-energy tensor of the source fields will be 

= P(o)W+'5p(x,t), 
T^' = (P(o)(t)+P(o)(i))5.0(x,t), (12) 

where 6 is the velocity potential of the fluid, 11 is the 
anisotropic stress tensor (shear) potential, Sp and Sp are 
small perturbations on top of the background homoge- 
neous density p(Q) and pressure P(q), respectively, and 
~ ^0 <C 1. We have used the Newtonian gauge because 
the gauge invariant scalar perturbations of the metric (<i> 
and ^) and gauge invariant perturbations in the stress 
energy tensor coincide with the perturbations explicitly 
written in this gauge. 

The general form of the equations for the perturbations 
in a metric theory with arbitrarily high derivatives is 



(13) 
(14) 



(15) 
I 

homogeneous equations. 

However, under a single assumption, it is possible to 
greatly reduce this freedom. The assumption is that 
the standard linearized equations of the perturbations 
of General Relativity are effectively recovered in the 
limit ^ H ^ for all the Fourier modes with 
i7+ ^ k > (subhorizon modes in the present time). 

Let us work out the implications of this assumption. 
In the equation for the perturbation of the energy density 
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in General relativity, 

8t:GSp= -6H^(j)-6H^+ ^A^: , (16) 

there is only a term proportional to Aip. However, for 
subhorizon modes {k ^ H) the dominant terms are those 
with the highest number of spatial derivatives. If we de- 
mand the terms with spatial derivatives of order greater 
than 2 not to spoil the behavior of these modes in the 
period in which GR is a good approximation, they must 
be negligible at least for the observable modes. 



For instance, they should be negligible for the modes 
responsible of the acoustic peaks of the CMB spectrum, 
and therefore for the modes which have entered the hori- 
zon after recombination (which have even lower k). It is 
possible that these terms are suppressed by inverse pow- 
ers of the UV scale becoming irrelevant for current 
tests of gravity, although they may be relevant for the 
physics of quantum fluctuations in the very early uni- 
verse. 

Therefore, at times when H ^ the equations can 
be approximated by 



SttGSp = 



B7rG(p(o) +P{o))di( 



m=0 
oo 

m—O 

oo 



bom _Q^^^O^^g. 



JJm-2 t ^ ^ JJm-2 t ^ ]-[m g2 i 



Corn ™ 



^TTG{-5p5] + d^djIl) = 5]Y^ 



m=0 

oo 

m=0 



"-Om 



bl,n A 



(17) 
(18) 



(19) 



where general covariance has been used to exclude the 
terms with coefficients eim, fim which are not compatible 



with (171, (181 



We must also take into account that the modes which 
are responsible of the acoustic peaks of the CMB un- 
dergo a phase in which they oscillate as sound waves, 
i.e.: dt(f> ~ kef). In order to explain the acoustic peaks of 
the spectrum of the CMB it is required that the acoustic 
oscillations of the modes of the gravitational potentials, 
0k and V'k, which lead to the acoustic peaks of the CMB 
spectrum, have a frequency ~ k. When radiation domi- 



nates and we consider modes well inside the horizon, two 
of the solutions of the system of differential equations of 
arbitrary order have this property ^. 

However, in order for nondecaying superhorizon modes 
to evolve into oscillating subhorizon modes with fre- 
quency ~ fc, a fine-tuning of coefficients is required unless 
the equations are of order two in time derivatives. Then 
the acoustic peaks in the CMB spectrum are reproduced 
when terms with more than two derivatives in the equa- 
tions for the perturbations can be neglected. Then 



A 



Oio- 



A 



bio—^ , 



m—O 



(20) 
(21) 



8^TG{-5pS) + ^^^JI^) = 



m=0 
2 

m=0 



a 



(22) 



r 



^ As we will see below, if the term with the highest time derivative is of order D, = -eoD as a consequence of general co- 
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and the arbitrariness in the evolution equations for the 
scalar perturbations has been reduced to twenty unde- 
termined dimensionless coefficients at this level. 



The equations (20 1,(21 1,(|22| are in principle valid for 
any perturbation mode {0k,V'k} as long as fc ^ H+. In 
particular, it must be valid for the mode k = 0, which 
corresponds formally to a perturbation with no spatial 
dependence. In practice, we will be able to neglect the 
spatial dependence of perturbation whose spatial depen- 
dence is sufficiently smooth, i.e.: its wavenumber k is 
sufficiently low. In GR it suffices for a mode to be super- 
horizon k ^ H m order to neglect its spatial dependence 
in a first approximation. 

If we consider the FRW metric perturbed by one of 
these modes we have, 



(1 + 2(j){t))dt^ - a^{t){l - 2ilj{t))dy? 



(23) 



>From now on the dot will represent derivative with 
respect to the cosmic time. By means of the invari- 
ance under time reparameterizations we can introduce a 
new time variable dt' = (1 -I- ip{t))dt, a new scale fac- 
tor a'{t') — a{t){l — 'ip{t)) and a new energy density 
p'(<') = po(t) + 5p[t) leading us back to the ACM in 
a homogeneous background. The Hubble rate H — ^ 
and its time derivatives iJ^") will change as 



dt'" 



da' 



a'(t') dt' 



n - 
m 



(24) 

Noticeably, this procedure must be applied to the gener- 
alized Friedman equation before one of its homogeneous 
solutions is used to write the time variable t as a function 
of H. Therefore we should start by considering a gener- 
alized first Friedman equation in primed coordinates 



8^Gp'(t') = V{H'{t'),H'{t'), i?(")'(i'), .. 



(25) 



and the corresponding equation for the pressure will be 
derived using the continuity equation. The linearized 
equations for the sufficiently smooth scalar perturbations 



in the Newtonian gauge ( 23 1 can then be derived with the 



help of (24 1. The result is 



87rG(5p(t) 



Y^g\^{H{t))5H^-\t), (26) 

^ oo oo 



H 



n— m— 

oo 



^ff|„(i/)(^. ^SH 



n=0 



n=0 



where 



9\niH) = 
g\nm{H) = 



df 



H("i=H('^)it{H)) 



/i-(n)=H(")(t(_ff)) 



(27) 

(28) 
, (29) 



and the time dependence of H is derived from the homo- 
geneous evolution This procedure fixes exactly the 



coefficients aom,6om,eor„ and /q™ in (17l,(18l,(19l, i.e. 



all the terms with no spatial derivatives in the equations 
for the scalar perturbations, as functions of f{H,H, ...) 
and its partial derivatives. 

We can now count the number of functional degrees 
of freedom in the superhorizon modes of the linearized 
theory. Let us assume that there is a maximum num- 
ber D of time derivatives in the equations for the per- 
turbations ^. The number of derivatives, with the use 



of (24 1, (26 1, (27 1, fixes / to be a function of at most 
When particularized to a solution of the homo- 
geneous Friedman equations, / and its first and second 
partial derivatives become functions of the Hubble rate 
H: g{H), g\n{H) and g\nm{H) respectively. That makes 
l-|-(i:)-l)-|-D(i:»-l)/2 = i:)(D-Hl)/2 functional degrees 
of freedom. The first derivative of g{H) can be written 
in terms of the g^n{H), and the first derivative of the 
later can be also written in terms of the g^nm{H). That 
makes 1 -f (_D — 1) conditions, so the result depends on 
D{D — 1) /2 independent functions of the Hubble param- 
eter H in the terms without any spatial derivatives. 

Notice however that if we restrict the number of deriva- 
tives appearing in the equations for the perturbations to 
two as in Eqs. ([20) , ([2lJ , (|22) , then Eqs. (|24) , ([26) , @ tell 
us that we must consider just functions f{H, H,H, ...) — 
f{H) = g{H), at least as an approximation at times 
H <C and modes k <C Therefore the number of 
functional degrees of freedom in the terms with no spa- 
tial derivatives is just one: the homogeneous evolution 
f{H) = g{H). 



variance. Well inside the horizon these are the dominant terms 
appearing in the wave equation resulting when calculating the 
adiabatic perturbations in the radiation domination epoch. 



^ We are forgetting here the restriction on the number of deriva- 
tives required in order to reproduce the acoustic peaks of CMB. 
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Our assumption might be relaxed and we could impose 
that general relativity should be valid for all the modes 
< k < kobs that have been observed in the spectrum 
of CMB and matter perturbations. This could lead to 
new terms in the equations for the perturbations (sup- 
pressed not necessarily by the UV scale H+) which have 
been negligible for the observed modes but that could 
lead to ultraviolet deviations from the spectrum derived 
in the general relativistic cosmology which have not yet 
been observed. However, terms with more than two spa- 
tial derivatives will be very tightly constrained by solar 
system experiments. 

In this article we will restrict ourselves to the system of 



us now derive the linearized equations for the rest of the 
modes under this assumption. In the (t',x) coordinate 
system, the metric is Friedman Robertson Walker, and 
we know that for this metric, we can use equation (|3]|. 
Thus, 



with 



(30) 



(31) 



second order differential equations (201, (21 1, (22 I. Let 



at linear order in perturbations. Therefore, we can de- 
duce from (|3| that in the (i, x) coordinate system 



1 

87rG(/9(o) -Hp(o))i9»6' = 5, ^ 



A A 



H 



8nG{-Sp6] + d,djU) = 6} 



H 



/\H){H<t> + i^) 



g'{H) (3i/0 + 3V;--^V'+§^ 



(32) 
(33) 



4>+ -j^i> 



(34) 



and the twelve coefficients of terms with no spatial deriva- 
tives are fixed by the function g{H) which defines the 
homogeneous cosmological model. 

The last requirement comes again from general covari- 
ance. If the tensor T^^ comes from the variation of a 
certain matter action Sm with respect to the inverse of 
the metric g'^", and Sm is a scalar under general coor- 
dinate transformations, then must be a 1-covariant, 
1-contravariant divergenceless tensor, i.e. V ^Tjf = 0. As 
the stress-energy tensor is proportional to G{^, the later 
must also be divergenceless. In the linearized approxi- 
mation = G|q-|j^ -I- SC^ and the Christoffel symbols 



version of this requirement as 



SV^jj , and we can write the linearized 



0. 



(35) 



Let us study this condition in order to see if we can 
further limit the number of independent coefficients of 
the equations for scalar perturbations. For ji — i the term 
(5G^|g will give at most a term proportional to doi^/'ii), 

which will not be present in other terms except SG^^^ 



This term will give at most a term g'ip^i/H which fixes 
the exact value of dgi for all H ^ . For the same 
reason, the term cqi^i^ in SG^^^ can not be canceled and 

the term cooH(j)\i can only be canceled by the term g'4)\i 
in SG^y, which means that cqi — and fixes the exact 
value of Coo for all H ^ The value of doo can be 

fixed in terms of doi ^ind the terms proportional to ^ in 



(34 1. This fixes (331 completely in terms of the function 
describing the homogeneous evolution, g{H): 



9'{H ) 
H 



(36) 



In the previous equation and from now on we will refer 
directly to the gauge invariant counterparts of the vari- 
ables in the Newtonian gauge, which we will denote by 
<&, ^, 5p and 5p (11 is already gauge invariant). Now we 
can use the divergenceless condition for /i = 0. The term 
(5Gq|q will give at most a term proportional to aioA0 and 

a term proportional to bio^V', which will not be present 
in other terms except <5Gg|^ . This term will give at most a 

term proportional to Ac/) and a term proportional to Ai/j. 
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Therefore aio = and the value of 610 is set completely 
in terms of g{H): 



8nG6p = -3g'{H){H^ + ii) 



9'{H) 



(37) 



The complete knowledge of (371 and ( 36 1 when H ^ 



gives a complete knowledge of (34 1 in terms of g{H) 
when H H+: 



i^TG{5p5^J+^^^J\l) 



H ^ 



S^, {g'iH) 

+ ^ {S,,A 1 d^d,) (^(<i> - VP) + 4,{g'{H) - Hg"{H))^ 



(38) 



and we have finally a set of equations for the scalar per- 



turbations (371, (36 1, (38 1 which are determined by the 



function g{H) which defined the ACM in the homoge- 
neous approximation. 

Had we relaxed our assumption, we could work with 



the system of equations (13l,(14l,(15l restricted to a 



maximum number of derivatives D as an approximation. 
Together with the condition G(^.^ = 0, the system of 
equations for the perturbations define a set of coupled 
differential equations for the coefficients of the terms with 
at least one spatial derivative. The set of equations com- 
ing from G^.^ = define the terms proportional to emn 
and fran in (15 I as a function of the terms in (14 1. The 



set of equations coming from Gq.^ define the rest of the 
terms in (151 as a function of the terms in (14 1 and (13 1. 



Therefore, the only freedom, for what concerns linearized 
perturbations, is that of choosing the set of functions 
{a„m, 6„m, c„m, d„m}, with some of them fixed by the 
homogeneous dynamics ([ij. An interesting case is found 
if D = 4 is imposed. This includes a description of f{R)- 
theories jl^, bigravity theories [20], and other of the most 
studied modified gravity theories. We will study further 
this case in the subsection devoted to the comparison of 
the model with /(i?) theories. 

The rhs of the equations derived are the equivalent to 
the components of the linearized Einstein tensor derived 
in a general covariant theory whose Friedman equation 
is exactly (|3}. Aside from General Relativity with or 
without a cosmological constant, i.e. for q;_ ^ 1,0, it 
will be necessary to build an action depending on arbi- 
trarily high derivatives of the metric in order to derive a 
theory such that both the equations in the homogeneous 
approximation and the linearized equations for the scalar 
perturbations are of finite order. 

Noticeably, the nonlinear equations for the perturba- 
tions will include derivatives of arbitrary order of the 
metric perturbations, with the exceptional case of Gen- 
eral Relativity with a cosmological constant. This makes 
problematic the consistency of the ACM beyond the lin- 
ear approximation. This issue will be further studied in 
a future work. 



To summarize, what we have shown in this section is 
that the behavior of scalar perturbations in a general co- 
variant theory is intimately connected to the background 
evolution. Except in the very early universe, the lin- 
earized equations for the scalar perturbations are deter- 
mined by the equations in the homogeneous approxima- 
tion if one assumes that there are no terms with more 
than two spatial derivatives. This assumption could be 
relaxed in order to include more general theories. 

There will be a subset of general covariant theories 
with a background evolution given by (|3} that verify 
these equations for the perturbations. We may be able to 
distinguish among these at the linearized level by means 
of the vector and tensor perturbations. 



HYDRODYNAMICAL PERTURBATIONS 

Vector perturbations represent rotational flows which 
decay very quickly in the General Relativistic theory. As 
we expect a small modification of the behavior of pertur- 
bations just in the vicinity of the lower bound H_, we 
will assume that the vector perturbations have decayed 
to negligible values when the scale begins to play a 
role and therefore they will be ignored. 

Tensor perturbations correspond to gravitational 
waves, which at the present have been not observed, their 
effect being far beyond the resolution of current observa- 
tions. 

Present measurements restrict their attention to per- 
turbations in the photon sector (CMB) and the matter 
sector (matter power spectrum, both dark and baryonic). 
These observations can be computed taking into account 
only the effect of scalar perturbations. Therefore, we 
can use the result of the previous section to study the 
deviations from ACDM in the spectrum of CMB. A com- 
prehensive study of the matter power spectrum predicted 
by the ACM would require a detailed knowledge of the 
nonlinear regime, which we lack at present. 

Let us start by deriving the behavior of scalar pertur- 
bations in a universe whose gravitation is described by 
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the ACM and filled with a perfect fluid (11 = 0) in ([12 I 



Therefore, (38 1 with j fixes $ as a function of ^ ^ 



H 
2IP 



1 - 



Hg"{Hy 
g'{H) , 



(39) 



Then, this expression can be used to substitute <I> in the 
remaining three differential equations for 5p, 5p and 0. 
Given the equation of state of the fiuid p = p{p, S), with 
S the entropy density, the perturbation of the pressure 
can be written as 



Sp = c^Sp + tSS , 



(40) 



where = (dp/dp) s is the square of the speed of sound 
in the fiuid and t = {dp/dS)p. Substituting the pressure 
and the energy density perturbations for the correspond- 



ing functions of ^ and its derivatives into (401, we arrive 



at the equation for the entropy density perturbation SS. 
The perturbations we are interested in are adiabatic, i.e.: 
dS = and therefore, the equation for the entropy per- 
turbations turns into an equation of motion of the adia- 
batic perturbations of the metric. If we want to consider 
the effect of the lower bound of the Hubble rate, on 
the evolution of perturbations, we must consider a matter 
dominated universe (c^ — 0), 



34/ - J,* 



H * 



0. 



(41) 



This equation is a second order linear differential equa- 
tion with non-constant coefficients. Thus it is useful to 
work with the Fourier components of the metric per- 
turbation ^fc. The equation of motion for the Fourier 
components is just a second order linear ODE with non- 
constant coefficients, which will be analytically solvable 
just for some simple choices of g{H). However, in general, 
it will be mandatory to perform a numerical analysis. 

THE CMB SPECTRUM IN ACM VERSUS ACDM 

The main observation that can be confronted with the 
predictions of the theory of cosmological perturbations at 
the linearized level is the spectrum of temperature fiuc- 
tuations of the Cosmic Microwave Background (CMB), 
from which WMAP has recorded accurate measurements 
for five years |21j . 

The temperature fiuctuations ^ are connected to 
the metric perturbations via the Sachs- Wolfe effect |22) . 
which states that, along the geodesic of a light ray 
^ = r(l + <f> + vp), characterized by the unit three- 
vector P, the temperature fiuctuations evolve according 
to 



Neglecting a local monopole and dipole contribution, 
taking recombination to be instantaneous at a certain 
time and assuming that the reionization optical depth is 
negligible, the present temperature fiuctuation of a dis- 
tribution of photons coming from a given direction of the 
sky r can be related to the temperature perturbation 
and metric perturbation ^ in the last scattering surface 
plus a line integral along the geodesic of the photons of 
the derivatives of <I> and "if (Integrated Sachs Wolfe effect, 
ISW). 

These relations are purely kinematical and remain un- 
changed in the ACM. In order to take into account all the 
effects involved in the calculation of the CMB spectrum 
the code of CAMB [23] has been adapted to the ACM. 
(see the appendix for details). 

We have compared the late time evolution of the per- 
turbations described by the ACM ^ for several values 
of a_ (taking into account that a_ = 1 corresponds to 
ACDM) for constant Hq and fi^ defined in Eq. (Il0|. We 



have assumed a Harrison-Zel'dovich scale invariant spec- 
trum of scalar perturbations as initial condition (n^ = 0) 
in the region in which the universe is radiation dominated 
but H ^ . Thus the study of the effect of the scale 
i7+ is postponed. We have taken Hq = 72 ± 8 km/sMpc 
from the results of the Hubble Key Project |2l] and 
rim = 0.26 ± 4 from our previous analysis of the evo- 
lution at the homogeneous level. The effect of ACM as 
compared to ACDM is twofold (FIG. 5). 



Tenperature Perturbations On = 0^36, H = 72 




FIG. 5: Spectrum of temperature fluctuations of the Cosmic 
Microwave Background for f2„i = 0.26, Hq = 72km/ sMpc 
and varying a_ = 0.3 (red), 0.6, 1.0, 1.5, 2.0, 3.0, 4.0 (blue) . 
The lowest values of a_ show extreme ISW effects and there- 
fore will be ruled out also by the CMB spectrum. Unfortu- 
nately, the cosmic variance masks the effect of 1.0 < a- < 4.0 
at large angular scales and therefore we can obtain little in- 
formation about the ACM from the late ISW effect. 



There is a shift in the peak positions due to the dif- 
ference in the distance to the last scattering surface for 
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different values of a-. This shift increases with the peak 
number. In particular, more precise measurements of the 
position of the third peak could be used to estimate the 
value of a_. 

There is also an increase of the lower multipoles due 
to the late Integrated Sachs- Wolfe effect, which is par- 
ticularly extreme for a_ ^ 1.0 (which corresponds to 
H- ~ Hq for constant fim)- This large deviation is 
clearly not present in the experimental data. The much 
smaller deviation for a_ > 1.0 is within the error bands 
due to the cosmic variance, and therefore the value of the 
lower multipoles cannot be used to exclude values of a- 
greater than but of the order of one. 



COMPARISON WITH THE TREATMENT OF 
PERTURBATIONS IN OTHER MODELS 

Perturbations in f(R) theories 

In our previous work jl3] we found that given an ex- 
pansion history parameterized by a modified Friedman 
equation in a universe filled by a given component, it 



was always possible to find a biparametric family of f{R) 
theories (see Ref. [27] for a review) which had the same 
homogeneous evolution as a solution of their equations of 
motion. We wonder now if these theories have also the 
same linearized equations for the scalar perturbations of 
the metric. 

An f{R) theory is defined by its action and therefore 
the equations for the perturbations are uniquely deter- 
mined. The action is given by 



S- 



1 



(43) 



where R is the curvature scalar and f{R) is an arbitrary 
function of this scalar. The equations governing the evo- 
lution of the perturbations will be derived from the Ein- 
stein equations of the f{R) theory, 

f{R)R^ - lf{R)6i^ + {Si^D - V^V,)/'(i?) = SnGTi: . 

(44) 

At the homogeneous level, it is possible to write the 
Friedman equations of the f{R) theory. 



8^Gp(o) = -HH + H^)f'iR^o)) - 2/(^(0)) - 3H.r(i?(o))^(o) , 
-87rGp(o) = -{H + 3H^)f'{R^o)) ~ ^/(i?(o)) + /('H^(o))^(o) + /"(^(o))i?(o) , 



(45) 
(46) 



where 



i?(o)(t) = -6(ij + 2ij2) 



(47) 



is the curvature scalar of the Friedman-Robertson- Walker 
(FRW) metric (zeroth order in perturbations). The so- 
lution of the system of equations (45 1,(46 1 together with 



the equation of state of the dominant component of the 
stress-energy tensor, defines i?(Q) and its derivatives as 
functions of the Hubble parameter H [13j. 



If we expand (44 1 in powers of the scalar perturbations 



of the FRW metric (111, the first order term gives the 
linearized equations for the perturbations. This will be 
a system of differential equations of fourth order, but it 
will be possible to turn it into a system of differential 
equations of second order, as we will see below. 

If the stress-energy tensor has no shear, the Einstein 
equation for 11 = i =/= i> = j reduces to 



f"{R(„))SR + f'{Ro){'P-^)^0, (48) 
where 6R is the perturbation of the curvature scalar 



SR = 12{H + 2ij2)$ + 6ij4 + 24Hi' 
+6<i'+ ^A($-2^')). 



(49) 



Eq. (|48| is a second order differential equation for the 
scalar perturbations except in the case of General Rela- 
tivity (/" = 0) where it becomes an algebraic equation 
($ — v[/). Eq. (48 1 can be used to turn the remaining 



components of the equations of f{R) theories (44 1 into 
second order equations. 

In the case of adiabatic perturbations in a matter dom- 
inated epoch one has 
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-SttG^P = /'(i?(o))(a-2A(4' - 2*) + ^ + 677* + + 2(77 + Si?^)^ - l5R) 
-d? [/'(i?(o))(<i> - - 2Hdt [/'(i?(o))(<i> - *)] = . 



(50) 



If we compare the equations governing the behavior of the other equations of the system (44 1. The equation for 



adiabatic perturbations in a /(i?) theory (48 1, (50 1 with 
those coming from the ACM, (391,(41 1, we find that they 



describe very different behaviors. In one case we have a 
system of two coupled second order differential equations 
for $ and \E'. In the ACM case we got a single second 
order differential equation and an algebraic equation be- 
tween the two gravitational potentials. For a given f{R) 
theory with an associated g{H) homogeneous behavior, 
we find that the behavior of linearized perturbations dif- 
fers from the one defined by the linearized perturbations 
in the ACM. This means that f{R) theories break in 
general our assumption in the third section (there are 
terms with more than two derivatives in the equations 
for the perturbations, and therefore the behavior of per- 
turbations in General Relativity is not recovered when 
i7+ > > i7_ for all modes) . 

Let us see this in a simple example. We will choose 
the easiest biparametric family of f{R) theories: the one 
that, under matter domination, gives the same back- 
ground evolution as General Relativity without a cos- 
mological constant |13] . 



f{R) = i? + ci|i?| A(7-v^) ^ ^^i^p 



r (7+773) 



(51) 



These kind of models were introduced in Ref. |28) . 

Let us first check if the solutions of the equations for 
the perturbations in General Relativity in the presence 
of pressureless matter. 



the pressure perturbations then gives 



-/'(^(o))[^- 
-27i/3/"(^(o))^- 



- 4i/$ + 3H^<i>] 
3$) = 8ttG6p ^ 



(54) 

It is obvious that the i ^ j equation of the cosmologic 
perturbations in General Relativity, 



$-^' = 



(55) 



is not recovered from (481 in the H 3> limit for all 
k ^ In fact the deviation would be significant for 

modes with k ^ Ha ( jp- 1 Ha. Therefore, 

the assumption we have made in order to derive the equa- 
tions for the perturbations in the ACM is broken. 

It is a known problem that f{R) theories are unable to 
pass cosmological and astrophysical tests involving per- 
turbations [29], unless the function involved is properly 
fine-tuned. In particular, if the conformal equivalence be- 
tween f{R) theories and scalar-tensor theories is used, it 
is necessary that the effective mass acquired by the new 
scalar degree of freedom is unnaturally large |30] . Some 
f{R) theories that pass cosmological and solar system 
tests have been proposed |31l l32l |33] . However, it is also 
subject of debate if the conformal equivalence can be used 
in order to extract physical predictions from this models, 
especially predictions which involve perturbations |34] . 
A mathematically rigorous treatment of perturbations in 
/(i?)-gravity can be found in Ref. |35] . 



^ — ^ = const , 
= $ cx t-^/^ , 



(52) 



Quintessence Models 



are solutions of the equations of motion in the case of the 
biparametric family of f{R) theories (51 1. The easiest is 
to verify if Eq. (44 1 with jj. ^ i ^ v ^ 



2/"(i?(o))(3i?^* + 3iJ$ + 12H^ + 3\^ + ^A($ - 2^-)) 
+/'(i?(o))('i>-*) = 0, 

(53) 

where R(q){H) = —3H^ is the value of the homogeneous 
curvature scalar as a function of H in this family of the- 



ories, is also fulfilled by the solutions (52 I. The result is 
obviously not. 

The second order differential equations for the pertur- 



bations in the f{R) theories ( 51 1 in the presence of matter 



can be found by substituting (53 1 and its derivatives into 



Another class of models used to describe the accel- 
eration of the universe are those in which a so called 
quintessence field, typically of scalar type (p, is added as 
a component of the universe [36] . The only functional de- 
gree of freedom in most models is just the scalar potential 
V{cp), which can be tuned to fit the homogeneous expan- 
sion of the universe. In our previous work jl3) we found 
the correspondence between a given homogeneous evo- 
lution parameterized by g{H) and the potential for the 
quintessence which drives under General Relativity this 
evolution. We now wonder if the behavior of perturba- 
tions in quintessence models also resembles the behavior 
in ACM. 

The field action is given by 
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(56) 



Let us assume that the field can be spHt in two compo- 
nents: one which is only time dependent and which drives 
the homogeneous evolution of the universe, (/3(o)(t), and 
a small perturbation which is inhomogeneous, 5(p(a;,i). 
We will consider linearized perturbations of the metric, 
the field, and the other components of the universe. The 
resulting linearized stress-energy tensor of the field is 



ip(o)5Lp - $y(i^(o)) + ]^V'{Lpi^Q))5ip , (57) 
(p(^o)dtS(p, (58) 

-fW(^(o)) - Jv'(^(o))'5^ . (59) 



By virtue of the Einstein's equations, it is always pos- 
sible to turn a modification in the Einstein tensor into a 
new component of the stress-energy momentum tensor of 
the sources of the gravitational field. We wonder if it is 
possible to account for the modification of the effective 
Einstein tensor described in the third section with a new 
component described by this quintessence field. However, 
it is not hard to see that the effective Einstein tensor that 
we are proposing has a modified i j component, while 
the i j component of the stress-energy tensor of the 
quintessence field (p is zero. Therefore, it is not possible 
to describe the evolution of perturbations in the ACM as 



driven by an effective scalar field component ( 56 1 . 



Dark Fluid Models 

A deformation of the gravitational physics can be also 
made equivalent to the addition of a non-standard fiuid 
component to the cosmic pie at the homogeneous level. 
The fiuid component used to explain the present acceler- 
ated expansion of the universe is typically taken to be a 
perfect fiuid with large negative pressure. 

One of the most popular parameterizations of this fiuid 
is the so called equation of state oj = p/p |37]- The sim- 
plest models are w = — 1, which is equivalent to a cosmo- 
logical constant, or a constant lu, but recently a possible 
time dependence of oj has been considered |38]. Need- 
less to say, these ad hoc parameterizations are subsets 
of the more general equation of state of a perfect fiuid, 
p — p{p, s), where s is the entropy density. It is also ques- 
tionable why must we restrict ourselves to perfect fiuids 
and not include a possible anisotropic stress tensor, pos- 
sibly depending on the energy and entropy densities |39) . 

As in the previous subsection, it is always possible to 
use the Einstein's equations to turn a modification in the 



gravitational physics into a new fiuid component of the 
universe. In the case of the equations derived in the third 
section, the equivalent Dark Energy Fluid would have the 
following properties: 



STvGdx 

SttGSpx 



m - g') 

m ~ g') 
H 



1 



A*-3(i?$ + 3') 



(60) 



(61) 



^{^('^-^^-sl^^^'-^^")*}'^^^) 

{2H - g'){3H^ + 3* - ^* + + ^ + 4*) 

+ ^{2~g"){H^ + ^)-&'KGMix. (63) 
ti 



The resulting fiuid is not an ideal fiuid (IIx = 0) or 
even a Newtonian fiuid {Jix (x 9x)- With the use of jH, 
(39 1 and (41 1, it will be possible to write 6px and 5px 
as a combination of Ox, and their derivatives. 



Spx = /i(P(o)x)Anx + /2(P(o)x)^x, (64) 
^Px = /3(P(o)x)nx - An^ + U{p(o)x)Ox ■ (65) 

These very non-standard properties show that, al- 
though it is formally possible to find a fiuid whose con- 
sequences mimic the ones of such a modification of the 
gravitational physics, this fiuid is very exotic. 

SUMMARY AND CONCLUSIONS 

In view of recent data, an updated comparison of cos- 
mological observations with a phenomenological model 
proposed in a recent work has been presented. 

An extension of this phenomenological model (ACM) 
beyond the homogeneous approximation has been intro- 
duced allowing us to describe the evolution of scalar per- 
turbations at the linear level. 

A comparison with the spectrum of thermal fiuctua- 
tions in CMB has been used to explore the possibility to 
determine the parameters of the ACM through its role in 
the evolution of scalar perturbations. The results of this 
comparison does not further restrict the parameters of 
the model, due to the masking of the associated late ISW 
effect by the cosmic variance. However, better measure- 
ments of the position of the third acoustic peak should 
improve the constraints significatively. 

It has been shown that the equivalence of different for- 
mulations of the accelerated expansion of the universe 
in the homogeneous approximation is lost when one con- 
siders inhomogeneities. In particular we have shown that 
the general structure of the evolution equations for scalar 
perturbations in the ACM differs from the structure of 
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the equations corresponding to modified f{R) theories of 
gravity, to quintessence models or to a dark fiuid with 
standard properties. 

The possibiUty of going beyond the Hnearized approxi- 
mation for the scalar perturbations and to consider vector 
and tensor perturbations will be the subject of a future 
work. 
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APPENDIX: COVARIANT PERTURBATION 
EQUATIONS FOR THE ACM 

The code of CAMB (23] makes use of the equations for 
the perturbations of the metric in the covariant approach. 
The quantities can be computed in a given "frame", la- 
beled by a 4- velocity u'^. In particular, CAMB uses the 
dark matter frame, in which the velocity of the dark mat- 
ter component is zero (the dark matter frame). Further- 
more, it parameterizes the time evolution with the con- 
formal time a(T)dr = dt. 

In order to apply CAMB to the ACM model it is nec- 



essary to identify frame invariant quantities, and then 
relate them to gauge invariant quantities [25]. The fol- 
lowing comoving frame quantities are used in the CAMB 
code: rj (the curvature perturbation), a (the shear 
scalar), z (the expansion rate perturbation), A (the ac- 
celeration, A — in the dark matter frame), <j> (the Weyl 
tensor perturbation), x'*-* (the energy density perturba- 
tion of the species i), g*^*^ (the heat fiux of the species i), 
and n^*' (the anisotropic stress of the species i). All of 
these quantities are defined in [2^. These variables are 
related to the gauge invariant variables via the following 
dictionary. 



-A+ = (67) 

= Sp, (68) 
k 

pq+{p + p)a = -{p + p)e, (69) 
a 

where prime denotes derivatives with respect to confor- 
mal time, except when acting on g{H) where it denotes 
a derivative with respect to the Hubble rate H, and 
T-C = a' /a ^ aH. On the other hand z = a+ -§^{1]' +2HA) 
and (f) — + ^)/2- The anisotropic stress 11 is already 
frame invariant. 

Written in the dark matter frame, the equations for 
the scalar perturbations read 



-nkg'{z' + nz) ^ B^GaE, H^il + 3ci'^^) - (H' - H)^il - ^) 



(70) 



The following combination of the constraint equations is also useful: 



k^(f> = - 



n 



(71) 



These equations are plugged into the Maple files pro- 
vided with the CAMB code and run to get the ISW effect. 
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